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SUMMARY 


A  limited  study  of  methods  of  deducing  generalized  air  loads 
acting  on  a  helicopter  rotor  blade  or  a  propeller  blade  from  bending- 
moment  measurements  is  presented.  A  new  set  of  basic  functions 
cal1  d  "orthonormal-moment  distributions"  is  proposed  as  the  best 
set  for  fitting  the  moment  data. 

The  blade -tip  deflections  computed  from  pressure  distributions 
measured  on  an  H-34  helicopter  rotor  in  steady  level  flight  at  112  knots 
are  compared  with  tip  deflections  derived  by  four  methods  based  on  use 
of  measured  bending-moment  distributions.  The  degree  of  agreement 
between  quantities  obtained  from  the  pressure  distributions  and  those 
obtained  from  the  moment  distributions  is  encouraging  but  not  satisfac¬ 
tory.  Neither  the  lift  nor  the  moment  distributions  were  measured  with 
sufficient  accuracy  for  the  purposes  of  the  present  study.  Deficiencies 
in  the  accuracy  of  the  measurements  and  in  the  accuracy  of  the  descrip¬ 
tion  of  the  blade  mass  ana  elastic  characteristics  obscure  any  meaning¬ 
ful  differences  between  the  various  methods  of  operating  on  the  moment 
data. 
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at  time  t  ,  pounds  0  ' 
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structural  area  moment  at  station  r  ,  inches 

equivalent  spring  constant  for  nih  mode  cu* J  m<f>*  dr) 
pounds  per  inch  °  ? 

local  lift,  pounds  per  inch 

measured  structural  moment  at  station  r  ,  inch-pounds 

maximum  mode  bending  moment  for  a  unit  generalized 
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natural  mode  moment  at  station  r  ,  inch-pounds 

orthogonal  moment  at  station  r  ,  inch-pounds 

mass  per  unit  spanwise  length  at  station  r' ,  stags" per  inch 

blade  tip  deflection  at  time  ,  t  ,,  inches 

blade  radius,  inches 

radial  coordinate  of  blade  measured  from  center  of' 
rotation,  inches 


time,  seconds 

part  of  aerodynamic  load  at  station  r  and  time  t  assumed 
independent  of  blade  response,  pounds  per  inch  (see  equation 
29) 

forward  flight  velocity,  inches  per  second 


IX 


(f! 

tj(r,  i) 


equivalent  mass  [=  J  m(r)  drj  ,  slugs 
bending  deflection  at  station  r  and  time  t 


inches 


i/h(^ 


bending  deflection  at  station  r  at  the  frequency  of  the 
fc^harmonic  of  rotational  speed  (fair)  =  yks(r)sin  kflt  +  t/kJr)cos 


e ■  error:  difference  between  measured  quantity  and 
approximating  series  at  station  L 

dfc  estimated  error  in  k^  harmonic  component  of  air  load 
measurement,  pounds  per  inch 

Ak  estimated  threshold  value  of  the  k^  harmonic  component 
of  the  air  load  measurement,  pounds  per  inch 


£  dummy  variable  of  integration 

2  4 

yO  density  of  air,  pounds  second  per  inch 

(j>n(r)  normalized  flapwise  bending  mode  shape  for  natural  mode, 
inches  per  inch  of  tip  deflection 


n  shaft  rotational  speed,  radians  per  second 


natural  frequency  for  mode  >7  ,  radians  per  second 
cO  frequency  of  oscillation,  radians  per  second 


azimuthal  angle,  radians 
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derivative  v/ith  respect  to  time:  -P 


df_ 

di 
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m,  n  bending  mode  number 
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INTRODUCTION 


Helicopter  rotors  and  yawed  propellers  are  subjected  to 
oscillatory  aerodynamic  loads.  These  arise  because  of  the  periodicity 
of  the  tangential  velocity  when  the  actuator  is  yawed  and  because  the 
wake,  being  unsymmetrical,  introduces  a  nonuniform  induced  velocity 
distribution  at  the  blading.  Theoretically  these  oscillating  loadings 
occur  at  all  frequencies  that  are  integral  multiples  of  the  rotational 
speed  in  the  steady-state,  yawed  flight  condition.  These  loads  at  the 
blading  feed  through  to  the  control  system  and  hub  so  that  not  only  the 
blades  but  the  linkages  in  the  control  system,  the  shafting,  and  the 
fuselage  as  well  are  all  subjected  continuously  to  periodic  loadings. 
Design  for  a  specified  fatigue  life  obviously  requires  some  knowledge 
of  these  oscillatory  loads. 

Few  direct  measurements  of  the  aerodynamic  forces  acting 
on  rotor  or  propeller  blades  exist,  and  only  recently  has  appreciable 
progress  been  made  in  the  development  of  theoretical  techniques  for 
predicting  the  higher -harmonic  air  loads.  Consequently,  blade 
designers  have  been  forced  to  extrapolate  past  experience  which,  for 
the  most  part,  has  been  based  on  interpretation  of  measured  bending- 
moment  distributions.  Rational  methods  for  interpreting  flapwise 
bending -moment  data  in  terms  cf  the  applied  loads  and  for  utilizing 
these  in  the  analyses  of  similar  configurations  were  developed.  The 
rationale  behind  these  methods  is  presented  in  References  1  and  2, 
and  these  expositions  are  highly  recommended  for  both  their  technical 
and  historical  importance. 

It  was  pointed  out  in  Reference  1  that,  if  perfect  bending - 
moment  measurements  were  made  and  if  the  mass  and  elastic 
characteristics  of  the  blade  were  known  exactly,  the  spanwise  aero¬ 
dynamic  load  distribution  could  be  obtained  from  measured  structural 
strains  to  within  the  approximations  of  simple  beam  theory.  The 
technique  would  involve  the  evaluation  of  differences  of  large  numbers 


and  the  double  differentiation  of  these  differences.  As  noted  in 
Reference  1,  the  accuracy  of  such  procedures  could  be  expected  to  be 
poor,  since  the  actual  data  are,  in  fact,  inexact.  It  was  believed  that  the 
combined  errors  resulting  from  measurement  inaccuracies,  approxima¬ 
tions  introduced  in  the  data  reduction, and  the  inherent  inaccuracy  of 
differentiating  the  experimental  data  would  make  the  calculation  of  the 
aerodynamic  loadings  meaningless. 

Since  this  direct  attack  on  the  problem  was  not  feasible, 
several  alternate  techniques  have  been  devised  witn  the  objective  of 
minimizing  at  least  those  errors  which  are  introduced  in  the  data 
reduction  procedures.  These  methods  generally  start  from  the 
assumption  that  the  azimuthal  variations  of  loads  and  bending  moments 
are  periodic.  The  bending -moment  data  are,  therefore,  harmonically 
analyzed  and  then  the  spanwise  variation  of  each  harmonic  is  assumed 
to  be  representable  by  the  sum  of  a  series  of  basic  functions.  These 
generic  functions  are  so  chosen  as  to  "most  likely"  fit  the  measured 
data  and  to  minimize  the  errors  in  the  subsequent  mathematical 
operations. 

Ail  of  the  analysis  procedures  which  have  previously  been 
proposed  used  one  or  the  other  of  two  basic  functions  —  either 
trigonometric  or  natural  mode  bending-moment  distributions.  Although, 
theoretically,  an  infinite  number  of  terms  is  required  to  define 
perfectly  the  moment  distributions,  only  a  few  seem  to  be  needed  in 
most  practical  cases  and,  in  fact,  the  number  is  limited  by  the  number 
of  data  points.  Variations  in  the  methods,  therefore,  generally  involve 
only  differences  in  the  number  of  terms  of  the  series  which  are 
retained  in  the  representation. 

The  method  which  uses  trigonometric  functions  as  the  basis 
is  treated  in  detail  in  Reference  3.  The  measured  moment  data  divided 
by  the  structural  rigidity,  El  ,  are  fit  with  a  truncated  Fourier  series 
whose  double  integral,  according  to  simple  beam  theory,  is  the  span- 
wise  deflection  curve..  This  curve  can  then  be  decomposed  into  a 


linear  combination  of  responses  in  the  natural-mode  shapes,  taking 
advantage  of  the  orthogonality  of  the  natural  modes.  CAL  has  utilized 
this  technique  to  calculate  blade  deflections  for  use  in  an  aerodynamic 
force  computation  (Reference  4). 

Approaches  which  use  the  natural-mode  bending -moment 
distributions  as  basic  functions  are  discussed  in  References  1,  2,  and 
5.  There  is  a  definite  spanwise  bending-moment  distribution  for  each 
natural  mode  which  can  be  obtained  from  beam  vibration  theory.  If 

XL 

the  transverse  deflection  of  the  blade,  ,  at  the  frequency  of  the  Jk 
harmonic  of  rotational  speed  is  represented  by  the  following  series  in 
the  natural  mode  shapes, 

oc 

y*.(r)  s  JZ  ?nJk  ’ 

n  =  i 

then  the  corresponding  bending-moment  distribution  is 

oo 

MTk(r)  =2T  mn(r)- 

n  •  i 

is  the  tip  deflection  due  to  bending  in  the  n*h  mode  at  the  ^  ^ 
harmonic  of  rotational  frequency.  The  <pns  and  7?lnS  are  considered 
to  be  known  quantities  determined  by  calculation  or  vibration  tests. 

A  particularly  simple  variation  of  this  approach  uses  only  a 
single  term  of  the  series  at  each  harmonic  frequency.  This  technique 
is  based  on  the  assumption  that  the  major  response  of  the  structure 
of  the  blade  at  each  frequency  is  in  the  natural  vibration  mode  nearest 
resonance  with  that  frequency.  All  other  modes  are  neglected,  and  the 
effects  of  interharmonic  coupling  for  the  particular  mode  are  neglected. 
This  single-mode  approximation  was  justified  on  the  basis  of 
experimental  observations  and  certain  physical  assumptions.  It  was 
further  argued  in  Reference  2  that  the  available  experimental  data 
contained  errors  of  the  same  order  as  those  made  in  neglecting 
intermodal  and  interharmonic  coupling.  This  approach  is  discussed 
in  detail  in  References  1  and  2  and  was  used  in  the  reduction  of  data 
reported  in  References  6  through  10. 
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Another  variation  which  utilizes  several  terms  of  the  series  in 
normal-mode  moment  distributions  rather  than  one,  as  in  the  preceding, 
is  utilized  in  Reference  5.  An  improvement  upon  the  single -mode 
approach  is,  in  this  case,  represented  by  consideration  of  the  first 
three  bending  modes  at  each  harmonic  frequency. 

None  of  these  techniques  have  proven  to  be  particularly 
satisfactory.  On  the  basis  of  physical  arguments,  the  normal-mode 
moment  distributions  constitute  a  more  satisfactory  set  of  basic 
functions  for  representing  the  total  moment  than  do  the  trigonometric 
functions.  The  neglect  of  the  intermodal  and  interharmonic  effects  in 
the  single-mode  approximation  seems  too  limited,  but  the  complexity 
of  the  calculations  increases  rapidly  as  the  number  of  modes  considered 
is  increased. 

Recent  events  provided  the  impetus  for  initiating  a  new  attack 
on  the  problem.  Successful  flight  of  helicopters  in  both  pure  and 
compound  modes  at  relatively  high  advance  ratios  showed  the  existence 
of  harmonic  excitations  at  flight  conditions  that  had  been  expected  to 
benefit  from  the  smoothing  introduced  by  a  large  freestream  velocity 
component  through  the  rotor  disc.  A  need  for  improved  prediction  of 
air  loads  was  obvious. 

A  basis  for  a  new  approach  resulted  from  the  observation  of 

& 

Dr.  H.  Daughaday  of  CAL  that  some  of  the  errors  introduced  in 
reducing  the  data  arose  frrm  the  lack  of  orthogonality  of  the  moment 
distributions  in  the  normf.1  modes.  His  suggestion  that  a  set  of 
orthogonal  basic  functions  might  be  generated  from  the  normal-mode 
moment  distributions  stimulated  the  development  of  the  method 
reported  here  based  on  "orthonormal-moment  distributions.  "  The 
desirability  of  using  orthogonal-moment  distributions  in  reducing  the 
data  is  based  on  theorems  of  linear  algebra  (see,  e.  g.  ,  Reference  11), 
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All  other  considerations  being  equal,  no  other  representation  will  be 
better  than  that  in  terms  of  orthonormal-moment  distributions. 

Perhaps  the  most  significant  event,  however,  was  the  in-flight 
simultaneous  measurements  of  pressure  distributions  and  structural 
moments  on  rotor  blades  in  a  series  of  programs  initiated  by  the  U.  S. 
Army  Aviation  Materiel  Laboratories  at  the  National  Aeronautics 
and  Space  Administration,  the  Bell  Helicopter  Corporation,  and  the 
Sikorsky  Aircraft  Corporation.  These  data  provide,  for  the  first  time, 
the  opportunity  for  comparing  results  obtained  from  moment  distribu- 
tions  with  those  obtained  from  the  measured  pressure  distributions. 

The  various  techniques  for  analyzing  the  bending -moment  data  might 
thereby  be  judged.  The  particular  quantities  compared  in  this  report 
are  the  tip  deflections  in  the  normal  modes  and  the  generalized  air 
loads  as  determined  from  the  measured  moment  and  pressure 
distributions. 

To  this  end,  the  derivation  of  the  modal  tip  deflections  from 
the  measured  lift  distributions  is  presented  in  the  first  section  and 
the  concept  of  a  "generalized  air  load"  is  introduced.  The  second 
section  contains  a  discussion  of  four  techniques  for  converting 
measured  bending  moments  to  modal  tip  deflections.  Values  of  tip 
deflections  obtained  by  the  various  bending-moment  data  reduction 
schemes  are  compared  with  the  corresponding  results  from  pressure 
distributions  in  the  third  section.  It  is  shown,  however,  that  the 
comparison  of  the  four  techniques  investigated  is  limited  by  the 
measurement  accuracy  of  the  available  data.  An  evaluation  of  the 
moment  measurer.  t  accuracy  required  to  perform  a  meaningful 
judgment  of  these  data  reduction  methods  is,  therefore,  presented 
in  the  fourth  section.  There  is  also  presented  in  this  section  a 
recapitulation  of  the  measurement  accuracy  estimate  given  in  the  data 
report  on  the  H-34  tests  (Reference  12), 
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It  should  be  recalled  at  this  point  that  the  objective  of  the 
studies  reported  in  References  1,  2,  and  5  was  the  development  of 
a  semi  empirical  method  for  estimating  blade  exciting  forces  for 
design  purposes.  Since  the  air-load  distributions  themselves  cannot 
be  determined  accurately,  the  concept  of  the  "generalized  force,  "  by 
means  of  which  the  effects  of  the  aerodynamic  loads  on  the  bending 
moments  might  be  approximated,  is  described  in  References  1,  2,  and 
5.  The  generalized  forces,  as  developed  in  these  references,  follow 
directly  from  the  modal  tip  deflections.  This  development  and  the 
important  assumptions  upon  which  it  is  based  are  reviewed  in  the 
fifth  section.  Generalized  forces,  as  defined,  do  not  include  the 
aerodynamic  forces  that  arise  from  bending  deformations  and  are 
discussed  in  relation  to  the  "generalized  air  loads"  which  include  the 
aerodynamic  damping,  intermodal,  and  interharmonic  forces. 
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1 .  THE  DETERMINATION  OF  MODAL  TIP  DEFLECTIONS 
AND  GENERALIZED  AIR  LOADS  FROM  MEASURED  LIFT 

DISTRIBUTIONS 


The  linearized  differential  equation  expressing  the  equilibrium 
of  transverse  forces  on  a  helicopter  rotor  blade  which  is  representable 
by  a  simple  beam  rotating  at  an  angular  velocity  fl  is 

2  /  a  -3.  aZ..  />r 


d 

3rz 


J (r-t-  §fl)-  (1) 


The  spanwise  distribution  of  aerodynamic  loading  is  represented  by  £ 
and  is  a  function  of  the  radial  position,  r  ,  the  time,  zf  ,  the  bla.de 
spanwise  bending  slope,  and  the  blade  bending  velocity.  The  bending 
deflection  of  the  blade,  y  ,  is  assumed  to  be  representable  by  a  double 
infinite  series  separable  in  r  and  t  which  reflects  the  periodicity  of  the 
loading  and  response  in  steady-state  forward  flight;  namely, 

OO  OP 

i/(r,6)  =21  ^n(r)X  ($n k.  coskflt  +$nk  sw  Mflt)  •  (2) 

r>  -  /  b=0 


<p„(rl  is  the  mode  shape  corresponding  to  the  blade  bending 
frequency,  u>n .  Each  of  the  bending-mode  deflection  shapes  satisfies 
equation  (1)  with  the  aerodynamic  loads  set  equal  to  zero,  since 
removal  of  the  aerodynamic  terms  reduces  the  equation  to  that  for  a 
vibrating  beam  in  a  vacuum.  Equation  (1)  then  becomes,  for  the  mode 
at  frequency  u)n> 


li- 

Sr* 


-.a 


f 

dr  \_dr 


r 


(3) 


Substitution  of  equation  (2)  into  equation  (1)  and  introduction  of  equation 
(3)  to  simplify  the  results  leads  to  the  following  equations  for  cosine 
and  sine  components: 
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o©  ~ 

Z  (cvn  -^CLZ)m<f>n  p„kc  =  Jkc 

(4) 

OO 

Z  (W -kZ£)Z)m<pn<{nk  =  ^ik 

o  =  /  *  s  s  J 

where  the  aerodynamic  loading  has  been  written  as 

J-Il-tkXr’%’W-)  cos  hm 

k -Q 

+EJ*s  (r’l£  >jj)  (5) 

h-0 

Since  the  natural  bending-mode  shapes,  <p ■■  ,  are  orthogonal,  the 
following  relationship  holds: 


J  m(r)  <pm(r)  <pn(r)c£r  =  O  fo 


r  m  ^  n  . 


Application  of  the  principle  of  virtual  work  and  the  orthogonality 
condition  permits  the  reduction  of  the  cosine  component,  for  example,  of 
equation  (4)  to 

■  Jc  JXr-7F'jf)  Wr)dr  -&-A-nkc  W 

where  (j^-j  -  j  m<f>n  is  the  equivalent  effective  mass  and  ^ nIc  is  the 

cosine  component  of  the  tip  deflection  of  the  nth  mode  at  the  kfh 

harmonic.  (r.A.n ^  is  defined  as  the  integral  over  the  span  of  the 

0 

product  of  the  nih  mode  shape  times  the  cosine  component  of  the  lift 
distribution,  ot  frequency  Lfl,  that  would  be  derived  from  measured 
pressure  distributions.  This  integral  will  be  called  the  "generalized 
air  load"  and  should  not  be  confused  with  the  integral  historically 
described  as  the  "generalised  force.  "  The  differences  in  these  two 
concepts  are  discussed  in  a  later  section  of  this  report.  For  the 
present,  it  is  sufficient  to  point  out  that  the  generalized  force  does  no* 
account  for  that  part  of  the  aerodynamic  loading  which  is  dependent  on 
the  blade  deformation. 
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From  equation  (6),  the  mortal  tip  deflections  are  proportional 
to  the  generalized  air  loads  in  accordance  with 


(7) 


2. 


THE  DETERMINATION  OF  MODAL  TIP  DEFLECTIONS 


FROM  MEASURED  BENDING -MOMENT  DISTRIBUTIONS 


Direct  measurement  of  blade  motions  by  photographic  means 
would  be  one  possible  method  for  obtaining  tip  deflections  from  which 
air  load  information  might  be  derived.  Generally,  however,  it  is 
simpler  to  measure  bending -moment  distributions  by  means  of  strain 
gages  bonded  to  the  blade  surface.  The  conversion  of  such  bending  - 
moment  measurements  to  modal  tip  deflections  is  the  subject 
of  this  section.  This  problem  is  essentially  one  of  curve  fitting,  and  it 
is  from  this  general  point  of  view  that  the  four  methods  described 
herein  are  related.  In  the  following  discussion  these  four  techniques 
are  called  (1)  Total-Deflection  Method,  (2)  Single-Mode  Method,  (3) 
Multi-Mode  Method,  and  (4)  Orthonormal-Moment  Method. 

The  general  technique  is  to  expand  the  measured  moment 
distributions  (or  functions  of  the  moment  distributions)  in  an  assumed 
series,  to  evaluate  the  coefficients,  and  to  deduce  deflections.  The 
mechanics  of  the  operations  are  sketched  in  the  discussion  below. 

Details  of  the  previously  employed  methods  can  be  found  in  References 
1,  2,  3,  and  5. 

The  four  techniques  to  be  discussed  differ  in  one  or  more  of 
the  following  characteristics:  (1)  the  quantity  fitted;  (2)  the  basis  used, 
i.  e. ,  the  functions  used  to  fit  the  quantity:  and  (3)  the  dimension  of  the 
basis,  i.  e.  ,  the  number  of  terms  retained  in  the  expansion.  As  pointed 
out  in  the  Introduction,  in  each  method  the  spanwise  variation  of  the 
bending  moment  is  represented  by  the  sum  of  a  series  of  basic  functions. 
The  methods  to  be  compared  utilize  as  basic  functions  either  trigono¬ 
metric  functions,  normal-mode  bending-moment  distributions,  or 
orthonormal  bending-moment  distributions.  The  first  of  these  is  used 
in  the  approach  which  will  be  called  the  "Total  Deflection"  method. 
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is  the  beam 


According  to  sxmple-beam  theory, 


£Z(iJ 


curvature  at  station  L  and  its  spanwise  variation  can  be  represented  by 

a  Fourier  sine  series  such  that  the  end  conditions  of  zero  moment  at 

the  root  and  tip  for  a  pinned-free  beam  (the  case  of  interest  here)  are 

satisfied.  This  series  is  integrated  twice  at  each  azimuthal  position 

(the  number  of  azimuthal  positions  being  governed  by  the  frequency 

range  of  interest),  and  two  constants  of  integration  appear  at  each 

dq 

position.  The  first  of  these  is  evaluated  by  equating  -r~  to  the 
measured  root  flapwise  slope,  and  the  second  constant  of  integration 
is  zero  to  satisfy  the  condition  of  zero  displacement,  y|  =0,  at  the 
root.  The  resulting  deflection  shapes,  are  decomposed  into 

modal  contributions  at  each  azimuthal  station,  if/p  ,  by  a  curve -fitting 
process;  that  is,  the  derived  deflection  shape  at  each  i//p  is  represented 
by  a  finite  series  in  which  the  natural-mode  deflection  shapes,  <pn (r)  , 

having  unit  tip  deflections  are  generic  functions.  The  deflection  shape 
is  represented  by 


d(r’  fp)  *p)  4>n  (  r) 


where  the  j>n(r)  are  known  from  calculations  or  vibration  tests  and  the 
coefficients,  ^.n(^p)  >  are  the  modal  tip  deflections.  The  mode  shapes 
weighted  by  the  square  root  of  the  mass  per  unit  span  are  orthogonal; 


i.  e. , 


rR 

J  m(r)  <pm(r)  4>n(r)  dr  =  0  ,  mfn 


so  that  equation  (8)  can  be  solved  for  the  tip  deflections: 


/  y(r,  4>p)  m(r)<pn(r)dr 

'  JL~p - T, - 

J  m(r)  4>n(r)dr 
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Under  the  assumption  of  periodicity  in  the  azimuthal  direction, 
expansion  of  yn  in  a  Fourier  series  in  yields  the  harmonic 
components  of  the  tip  deflections: 


(10) 


In  both  the  second  and  third  methods,  the  bending  moment 
data  are  first  harmonically  analyzed  with  respect  to  the  azimuthal 
variation.  The  spanwise  bending-moment  distribution  at  each 
harmonic  frequency,  kfl ,  is  tl  en  represented  by  the  sum  of  a  series 
of  normal-mode  bending-moment  distributions: 

oo 

fnttoj’-)  (11) 

n  =  f 


where  the  7?ln(r)  are  known  from  calculations  or  vibration  tests. 

In  the  approach  called  the  "Single-Mode"  method,  it  is 
assumed  that  all  the  bending  response  occurs  in  a  single  natural 
bending  mode  — that  estimated  to  be  near  resonance  at  the  frequency 
in  question.  Consequently,  the  measured  moment  distribution  is 
approximated  by  one  term  of  the  series  with  natural-mode  moments 
as  a  basis  [  equation  (11)].  In  effect,  the  spanwise  variation  of  the 
moment  at  each  harmonic  frequency  is  assumed  a  priori,  and  only  the 
amplitude  of  that  shape  remains  to  be  specified.  If  the  natural-mode 
moments  are  defined  ior  a  unit  tip  deflection,  then  the  coefficient,  to 
be  determined,  of  the  natural-mode  moment  distribution  is  the  tip 
deflection.  In  equation  form, 

MTk(r)  ^  m„Ar) 

or 
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Mtu  (r) 

Kc 


(r) 


(12) 


^T,  (r) 

a  /  5 

rn  k-s  ~  7T)  „<  (r) 

n  s 

The  subscript  n'  denotes  the  particular  natural  mode  whose  resonance 
is  nearest  the  harmonic  frequency,  k  fl  ,  being  considered. 


The  approach  called  the  "Multi-Mode"  method  also  proceeds 
from  an  expansion  of  the  moment  distribution  in  a  series  with  #2  (r) 
as  the  basis  and  coefficients  <?  .  In  this  case,  however,  several 

Ttlk 

terms  of  the  series  are  retained;  for  example, 

H 

*rt(r)  ink  *"„(?). 

n-  f 

The  dimension  A/ is  selected  such  that  the  natural  frequencies  of  the  N 

modes  completely  span  the  frequency  band  up  to  the  highest  harmonic 

of  shaft  rotational  speed  of  interest.  For  example,  if  the  highest 

harmonic  to  be  considered  is  the  tenth,  the  dimension  can  be  chosen 

as  four  since,  generally,  the  first  three  flapwise  bending  modes  have 

natural  frequencies  less  than  the  tenth  harmonic  and  the  fourth  mode 
frequency  is  not  far  above  it.  Of  course,  N  can  be  no  greater  than  the 

number  of  valid  measurements  of  moment  over  the  span. 


The  new  technique  developed  in  the  current  study  is  called  the 
"Orthonormal -Moment"  method.  This  approach  differs  from  the 
Single-Mode  and  Multi-Mode  methods  only  in  the  basic  function  which 
is  used  in  the  approximating  series  of  equation  (J  1).  In  all  other 
respects,  the  procedures  for  determining  the  modal  tip  deflections  are 
similar  in  these  three  methods.  A  new  orthogonal  set  of  basic 
functions,  ‘fon(r)  ,  is  derived  from  the  natural-mode  moment 
distributions. 
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The  ’desirability  of  using  orthogonal  moment  distributions  in 
fitting  the  measured  moment  data  is  based  on  the  theorems  of  linear 
algebra  (see,  e.  g.  ,  Reference  11),  Consider  the  measured  moment 
distribution  to  be  a  vector,  Mr  ,  and  the  normal-mode  moment 
distributions,  ,  to  be  linearly  independent  vectors  in  a  linear  mani¬ 
fold  belonging  to  the  same  space  as  the  measured  moment.  Then 
construct  an  orthonormal  set  of  vectors,  ,  which  span  the  same 
linear  manifold  as  the  to;  .  The  vector 


^TP(r)  =  2. 


r  1 


[  MT(Z) 


ht(r) 


(13) 


is  the  projection  of  MT  on  the  manifold  containing  the  ?&;  and  and  is, 
of  all  vectors  in  the  manifold,  the  nearest  to  Mj  . 


While  the  normal-mode  shapes,  <Pn  ,  form  an  orthogonal  set, 
the  corresponding  moment  distributions,  ???;  ,  in  general,  do  not,  since 
the  El  distribution  is  typically  not  constant  across  the  span  of  a  rotor 
blade.  It  is  possible,  however,  to  construct  an  orthonormal  set  of 
moments  from  the  natural-mode  moments  by  the  Gram-Schmidt 
process  (see,  e.  g.  ,  Reference  11).  A  vector  of  the  set  to  be  cc<  /erted 
(i.  e. ,  an  )  is  chosen  as  the  starting  point  and  the  first  orthonormal 
vector  is  constructed.  In  the  present  case,  the  moment  distribution,  \  , 
for  the  first  pinned -free  mode  is  a  convenient  starting  point,  and  the 
first  orthonormal  function  is 


7n,(r) 


7n,t(r) 


R. 


(14) 


The  second  orthonormal  function  is  constructed  from  the  second-mode 
moment  distribution,  7nz  ,  a  correction  term  to  make  ^perpendicular 
to  7nf ,  and  a  normalizing  divisor, 
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7rc?  (r)  = 


7nz(r)-  m^r)  £  7nt(Z)7nz(s)cL$ 
t If  *  [&,(*)- f  fn,  (4)  7nz(f)d?  dv 


(15) 


This  process  is  continued,  each  function  771  (r)  being  made  to  depend 
on  7nz(r)  and  to  be  orthogonal  to  all  prior  tH  .  That  is, 


i: 


7n;  ( r )  ?7i  (r)  oir  = 


(16) 


Now  the  measured  moment  distribution  can  be  written  as  a  series  in 
the  constructed  functions: 


Mr(r)  -JT  d-  7nL(r),  (17) 

6  -  / 

For  a  fixed,  finite  number  of  terms,  there  is  no  other  series  which 
would  constitute  a  better  fit  to  the  moment  measurements.  The 
amplitude  coefficients,  ,can  be  obtained  explicitly  from  the  measure¬ 
ments  as 


/?  _ 

>/>  (r)  77h-  (r)  dr 


(18) 


by  taking  advantage  of  the  orthonormality  of  the  %  functions. 


The  d-i  coefficients  can  be  related  directly  to  the  tip  deflections 
in  the  normal  modes.  Use  of  the  associated  matrix  method  (Reference 
13)  produces  a  distribution  of  moments  corresponding  to  a  unit  tip 
deflection  in  the  particular  mode  calculated.  Therefore,  the  measured 
moment  distribution  can  also  be  expanded  in  the  (nonorthogonal) 
natural -mode  moment  distributions 

co 

MT(r)  =  2;  Vlj(r)  (19) 

J- 1 

where  a,  is  the  tip  deflection  in  the  mode.  The  identity  of  the 
moment  in  the  two  coordinate  systems  [equations  (17)  and  (19)]  can 
be  written 


oo 

m/r)  cijMj(r). 


(20) 


Now  it  is  noted  that,  in  the  Gram-Schmidt  process,  the  7 U  functions  are 
linear  combinations  of  the  HI  functions,  so  they  can  be  expressed  as 

d7hi 


a7n» 


Mi 


so  that  equation  (20)  c  n  be  written 


(21) 


Rearrangement  of  the  right-hand  side  yields  the  following  convenient 
form: 


and,  by  equating  coefficients  of  , 


OO 


-  L 

‘’--j 


(22) 
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The  c£>;'s  are  evaluated  by  equation  (18)  and  the  - — i-  from  the 
* 

construction  of  the  orthonormal  functions  as  represented  by  equations 
such  as  (14)  and  (15).  The  latter  quantities,  ^  ..c  ,  can  be  evaluated 
as  accurately  as  desired, depending  on  the  number  of  stations  at  which 
the  normal-mode  bending  moments  are  evaluated  or  measured. 

In  all  but  the  Single-Mode  method,  there  is  a  choice  to  be 
made  of  the  number  of  terms  to  retain  of  a  particular  series 
representation.  As  pointed  out  previously,  the  number  can  never  be 
greater  than  is  consistent  with  the  number  of  valid  measurements.  If 
the  dimension  of  the  basis  selected  is  equal  to  the  number  of  measure¬ 
ment  points,  the  coefficients  in  the  approximating  function  can  be 
determined  by  requiring  that  the  approximating  function  pass  through 
the  measured  values.  If  the  dimension  is  less  than  the  number  of 
measurement  points,  the  coefficients  in  general  cannot  be  so  chosen. 

In  view  of  the  desire  to  minimize  the  number  of  terms  of  the  series 
retained,  it  is  this  latter  case  that  is  of  particular  practical  interest 
here.  The  difference  between  the  measurement  at  station  i  and  the 
approximating  function  at  station  l  is  defined  to  be  the  error,  &c  ;  i.  e, , 

ti  =  [^a)  -  hd] 

where 

F(i)  is  the  value  to  be  fitted  at  the  ^  station, 

Hi)  is  the  value  of  the  approximating  series  at  the  station. 
Generally,  some  criterion  related  to  the  error  is  used  to  permit 
evaluation  of  the  coefficients  of  the  base  vectors  in  the  approximating 
function  -P(i)  .  In  fact,  minimization  of  the  sum  of  the  squares  of  the 
errors  was  used  as  the  criterion  in  all  four  techniques  in  this  study. 

In  the  Total-Deflection  method,  the  least-square-error 
criterion  is  applied  to  the  fitting  of  the  generated  deflection  shape  by 
the  sum  of  a  few  normal-mode  shapes,  whereas  the  bending-moment 
measurements  themselves  are  approximated  in  the  other  three  method?. 
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In  the  Single-Mode  method,  the  error  is  defined  as 


~  ~  $n’k  7n 

and  minimization  of  the  sum  of  the  squares  of  these  errors  results  in 
the  following  equation  for  the  tip  deflection: 

p 

ZL  7nn/(c.) 

?«'*  «  - —  (23) 

Efau] 

0=1 

where  P  is  the  number  of  spanwise  measurement  stations. 

The  error  in  the  Multi-Mode  method  is  defined  as 

N 

~  ~£f  .fnt  771  n 

where  N  is  the  number  of  normal  modes  to  be  considered  and  is 
less  than  P .  Minimization  of  the  sum  of  the  squared  error  produces 
N  inhomogeneous  equations  from  which  the  tip  deflections  are 
determined. 


p  p  p 

" 

P 

Z  Kwj  Z7nt(i)7n^i)  *  *  * 

H 

ZMrZ)?n>a) 

0-1  0-1  0=1 

0  =  1 

p  p  z  p 

p 

•  •  *  Znza)mN(o) 

C r /  0-1  0-1 

m  *  « 

$Zk. 

Z  MrSl)mz  (L) 
0=1  ** 

•  •  • 

P  P  P  Z 

• 

p 

Z^Dm^u)  Z^K(t>  *  ‘  ’  Z  Lm»(L)] 

o=!  i  =  1  o  =  1 

%Nk. 
_  _ 

Zy^r,iN(i) 

This  operation  can  be  considerably  simplified  in  the 
Orthonormal-Moment  method.  As  in  the  Multi -Mode  method,  the 
error  is 

N 

»  MT  (i)-£cint(i)ihn(i) 

0  *■'  „  y 

n  b  i 

and  the  coefficients  are  determined  such  that  the  sum  of  the  squares 
of  these  errors  is  minimized.  The  set  of  equations  for  these 
coefficients  is  similar  to  that  in  equation  (24).  If  the  measurement 
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_  —  -  _ 


i 


stations  are  nearly  equally  spaced,  however,  the  following  simplifi¬ 
cation  can  be  introduced  as  a  consequence  of  the  orthogonality 

of  the  functions  T7l  : 

n 

P  _  .R  f  O  ;  m  f  n 

H  7rlm^7r)n(-L)  555  Co^sia.nt-J  Trirn(r)  ?nn(r)dr  =  < 

o-f  0  (  consta-n  t  >  m  =  r>  . 

The  equations  for  the  then  become  simply 

p 

II  MrkU)  m„U) 

*nt  ,  - -  (25) 

I  fJ*nY 

t -t 

The  coefficients  are  related  to  the  normal -mode  tip  deflections  by 
equation  (22).  The  above  factors  are  summarized  in  Table  I. 

TABLE  I 

SUMMARY  OF  CHARACTERISTICS  OF  METHODS  FOR  FITTING  MOMENT  DATA 


METHOD 

QUANTITY 

FITTED 

BASIC 

FUNCTION 

DIMENSION 

ERROR 

CRITERION 

"TOTAL 

DEFLECTION" 

Zp(r) 

DETERMINED  BY  NUMBER 

OF  MODES  HAVING  NATURAL 
FREQUENCIES  IN  FREQUENCY 
RANGE  OF  INTEREST. 

"SINGLE  MODE" 

A Ar  (D 

TWyjLr) 

1 

"MULTIPLE  MODE" 

"ORTHONORMAL 

MOMENTS" 

Mr  (i) 

t 

Tn  (r) 
n 

nn  (r) 

J 

DETERMINED  BY  NUMBER 

OF  MODES  HAVING  NATURAL 
FREQUENCIES  IN  FREQUENCY 
RANGE  OF  INTEREST. 

. 

kMIN IMUM 

2  s' 

Mr^(i)  =  Measured  moment  at  the  k?*'  harmonic  at  station  t  . 

7n,n(r)=-  Radial  distribution  of  natural-mode  moment  per  unit  tip  deflection. 
2>z-„(r)  =  Rad  ial  distribution  of  constructed  orthonormal  moments. 


#r  - 
El 


-  2 L  hj  sen 


.  (2,  -l)rr r 


;  where  P  ~  number  of  measurement  stations. 


I  JSC 
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It  is  appropriate  at  this  point  of  discussion  to  introduce  a  few 
comments  with  regard  to  consistent  da*a  processing.  The  least-square- 
error  criterion,  which  has  been  used,  weights  all  the  measurement 
errors  equally.  Practical  considerations  of  the  actual  measuring 
system,  however,  justify  the  intelligent  filtering  of  these  data  prior  to 
analysis.  For  example,  very  small  moment  amplitudes  should 
probably  be  discarded  on  the  basis  of  threshold  limitations  of  the 
measuring  system.  Moment  measurements  which  indicate  discontinu¬ 
ities  can  usually  also  be  discarded  as  being  physically  unrealistic. 

The  minimum  number  of  azimuthal  measurement  points  that  remains 
determines  the  highest  harmonic  that  could  possibly  be  considered;  the 
minimum  number  of  radial  stations  that  remains  determines  the  highest 
order  of  the  normal  or  orthonormal  mode  in  the  spanwise  representation. 
In  addition  to  the  threshold  considerations,  a  further  restriction  results 
from  the  limits  of  overall  accuracy  of  the  measurement  system.  Any 
mode  whose  maximum  amplitude  is  found  to  be  less  than  the  possible 
error  due  to  measurement  inaccuracy  should  properly  be  discarded. 

These  considerations  were  introduced  in  the  treatment  of  the 
data  used  in  the  next  two  sections. 
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3.  COMPARISON  OF  TIP  DEFLECTIONS  AND  GENERALIZED 
AIR  LOADS  DERIVED  FROM  MEASURED  LIFT  DISTRIBUTIONS 

AND  BENDING  MOMENTS 


It  was  pointed  out  in  Section  1  that  modal  tip  deflections  can 
be  obtained  from  the  generalized  air  loads  calculated  from  measured 
lift  distributions.  Methods  for  reducing  measured  bending-moment 
data  to  obtain  modal  tip  deflections  were  discussed  in  Section  2. 
Numerical  results  based  on  Data  Table  21  of  Reference  12  are 
presented  and  discussed  in  this  section. 

The  data  from  Table  21  of  Reference  12  which  were  used  in 
the  present  calculations  are  reproduced  in  Table  II.  The  aerodynamic 
load  distributions  presented  in  part  (b)  of  Table  II  were  used  to 
evaluate  the  generalized  air  loads  for  the  first  four  natural  bending 
modes.  The  results  are  presented  in  Table  III.  The  smaller  values 
indicated  in  Table  III  (that  is,  values  less  than  about  ten)  are  probably 
without  meaning  because  of  the  errors  in  the  measurement  and  the 
errors  introduced  in  the  integration.  Errors  in  the  integration  arise 
because  the  spanwise  distribution  of  measurement  stations  is  such  that 
curves  drawn  through  the  data  points  are  not  well  defined.  This  can  be 
seen  in  Figure  1  where  the  harmonic  air  load  data  are  presented.  The 
estimated  accuracy  of  the  values  indicated  in  Table  III  is  presented  in 
Section  4  as  part  of  an  effort  to  establish  corresponding  bending-moment 
accuracy  requirements. 

The  four  analysis  techniques  described  in  Section  2  were 
applied  to  the  bending -moment  measurement  data  presented  in  part 
(a)  of  Table  II. 

Table  IV  contains  a  listing  of  the  harmonics  of  modal  tip 
deflections  deduced  from  bending  moments  and  from  lift  distributions. 
The  deflections  corresponding  to  the  lift  distributions  are  used  as  the 
standard  in  the  following  discussion.  It  can  be  noted  immediately  that 
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TABLE  n 

FLIGHT  TEST  DATA  (REFERENCE  12) 


FLIGHT  18,  TRIM  LEVEL  FLIGHT  OUT-OF-GROUND  EFFECT  (V  =  112  KNOTS) 


(a)  Flapwise  bending  moment 
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TABLE  H 

COMPARISON  OF  DERIVED  HARMONIC  COMPONENTS  OF  MODAL  TIP  DEFLECTIONS 


the  results  based  on  any  of  the  four  methods  for  reducing  the  moment 
da^a  are  not  in  good  agreement  with  those  based  on  the  air  loads. 
Consequently,  a  detailed  inspection  of  these  results  is  not  warranted 
at  this  point.  Some  general  observations  can  be  made,  however. 

Tip  deflections,  rather  than  the  generalized  air  loads,  have 
been  tabulated,  since  the  physical  significances  of  their  relative 
magnitudes  are  more  easily  visualized.  According  to  the  results 
shown  in  Table  IV  based  on  the  air-load  measurements,  the  first¬ 
mode  displacements  dominate  at  the  first  three  harmonics,  and  first 
and  second  modes  dominate  at  the  fourth  and  fifth  harmonics. 

Deflection  amplitudes  are  indicated  to  be  small  (less  than  0.  1  inch) 
for  harmonics  above  the  fifth. 

Roughly,  the  tip  deflections  obtained  from  the  air-load 
analysis  and  the  bending-moment  analyses  follow  similar  trends. 

The  results  from  the  bending-moment  calculations  consistently 
show  larger  deflections  than  were  calculated  from  the  lift  distributions 
for  the  fourth  mode.  The  calculations  made  using  the  total  deflection 
technique  tend  to  indicate  larger  modal  deflections  than  any  of  the  other 
calculation  methods  for  modes  above  the  first  in  harmonics  above  the 
sixth. 


Although  tip  deflections  are  easily  visualized  they  are  not 

the  fundamental  quantities.  It  is  of  interest  to  convert  them  to  the 

corresponding  generalized  air  loads.  These  are  presented  in  Table  V. 

The  same  comparative  trends  exist  as  in  the  tip  deflection  data  (Table 

1 

IV):  but  the  effect  of  the  undamped  amplification  factor,  y—  ~z — TflTTiWY 

L  n  -  &  fl  )n 

is  to  make  the  relatively  small  higher -mode  tip  deflections  yield 
relatively  large  generalized  forces.  For  example,  using  the  tip 
deflections  derived  from  air  loads,  the  fourth-mode  ccsine  deflection 
is  less  than  1%  of  the  first-mode  cosine  deflection  at  the  first 
harmonic,  while  the  fourth-mode  cosine  generalized  air  load  is  12% 
of  the  first-mode  cosine  generalized  air  load  at  the  first  harmonic. 
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COMPARISON  OF  DERIVED  GENERALIZED  AIR  LOADS 


OA..k  CA.2L  GA.Jk  GA.  if.t 

Harmonic  _ 

Method  of  Order  (pound) _ (pound)  (pound)  (pound) 


Analysis 

t  -  7l 

Cos 

Sin 

Cos 

Sin 

Cos 

Sin 

Cos 

S,n 

Ai r  Load 

488.7 

-507.7 

-110.3 

141.6 

-148.  1 

139.7 

62.0 

14.4 

Total  Deflection 

529.4 

-478.4 

4.3 

-20.0 

70.4 

-28.5 

438.6 

-746.0 

Single  Mode 

1 
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-502.9 

- 

- 

- 

- 

- 

Multi  Mode 

446.5 

-513.2 

-31.3 

-32.0 

-134.2 

125.9 

-1 12.4 

•.£.7 

Ortho  Moments 

436.5 

-51  1.7 

-30.8 

-41.8 

-133.5 

127.0 

-114.4 

47.  ’ 

Ai r  Load 

-136.6 

139.0 

-92.6 

-IG7.3 

-40.6 

-57.2 

-73.9 

35.5 

Total  Deflection 

-143.2 

89.5 

-167.5 

-177.  1 

-183.  ! 

-23.3 

- 109'4.0 

-462.2 

Si ngl e  Mode 

2 

-90.3 

103.6 

- 

- 

- 

- 

- 

- 

Multi  Mode 

-102.7 

98.4 

-99.8 

-89,8 

42.5 

16.8 

265.8 

21.5 

Ortho  Moments 

-102.0 

98.9 

-101.5 

-91.5 

46.6 

19.9 

267.8 

22.5 

Ai  r  Load  . 

42.8 

-103.8 

-30.8 

10.2 

-15.6 

2.4 

-19.4 

-34.0 

Total  Defection 

43.7 

-67.2 

-77. S 

-33.7 

-  106.4 

-92.6 

-262.6 

-419.  1 

Single  Mode 

3 

21.7 

-71.4 

- 

- 

- 

- 

- 

- 

Multi  Mode 

19.0 

-71.5 

-47.2 

12.8 

-39.4 

21.9 

-4.7 

86.2 

0;  tho  Moments 

19.0 

-71.6 

-47.5 

12.9 

-39.  1 

22.2 

-4.7 

86.2 

Air  Load 

47.5 

-11.4 

-21.0 

-18.8 

-12.5 

14.8 

-6.4 

10.2 

Total  Deflection 

75.5 

-14.4 

-31.5 

-6.3 

-91.1 

66.7 

-180.6 

103.5 

Single  Mode 

1 

- 

- 

-27.7 

-7.3 

- 

- 

- 

- 

Hulti  Mode 

56.8 

9.7 

-23.4 

-3.6 

-6.8 

37.4 

-17.2 

47.2 

Ortho  Moments 

57.2 

10.9 

-23.3 

-3.9 

-7.3 

30.8 

-18.  i 

47.2 

Air  Load 

37.2 

4.8 

-3.  1 

8.2 

-11.6 

-1.  1 

4.7 

-5.6 

Total  Deflection 

14.5 

9.6 

-1.0 

1.3 

43.3 

35.4 

7.7 

116.3 

Single  Mode 

5 

- 

- 

1.  1 

1.9 

- 

- 

- 

- 

Hulti  Mode 

121.6 

18.9 

0.8 

1.9 

15.9 

-6.3 

-12.0 

2.6 

Ortho  Moments 

121.8 

19.3 

0.8 

1.9 

15.4 

-6.7 

-12.3 

2.6 

Ai r  Load 

-1.3 

37.3 

-12.4 

3.2 

3.7 

-21.8 

3.0 

-5.2 

Total  Deflection 

-26.  1 

49.  1 

-21.0 

21.5 

30.2 

-39.6 

94.2 

-138.5 

Single  Mods 

6 

- 

- 

12.4 

5.9 

- 

- 

- 

Hulti  Mode 

131.6 

6.3 

0.7 

5.9 

1.2 

2.9 

-16.6 

-7.  I 

Ortho  Moments 

131.5 

6.4 

0.7 

5.9 

1.2 

2.9 

-16.6 

-7.  1 

Ai r  Load 

13.0 

1.9 

-4.  1 

-2.3 

-1.6 

5.2 

9.4 

-0.5 

Total  Deflection 

1 17.  1 

38.4 

39.3 

24.  1 

-17.6 

-15.0 

-70.  1 

-30.8 

Single  Mode 

7 

- 

- 

- 

- 

0.4 

-8.3 

- 

- 

Multi  Mode 

150.9 

-74.  1 

2.8 

4.8 

2.4 

-i  1.2 

-10.0 

-8.6 

Ortho  Moments 

150.9 

-73.3 

2.9 

4.8 

2.4 

-11.2 

-10.7 

-8.6 

Ai r  Load 

-0.7 

10.8 

-1.7 

3.5 

-1.6 

-0.5 

-1.3 

-4.6 

Total  Deflection 

64.8 

-142.9 

21.1 

-71.6 

2.2 

-3.6 

-16.4 

121.4 

Single  Mode 

8 

- 

- 

- 

- 

0.9 

1.  1 

- 

- 

Multi  Mode 

125.7 

-154.9 

7.4 

5.  1 

SB 

1.8 

-13.2 

-30.2 

Ortho  Moments 

125.4 

-153.6 

7.6 

15.4 

HU 

1.8 

-13.2 

-30.2 

Ai r  Load 

4.6 

5.5 

0.7 

-5.2 

-2.4 

2.8 

-2.6 

-2.9 

Total  Deflection 

-139.6 

-59.4 

-116.8 

-18.7 

-31.3 

-3.7 

105.6 

40.3 

Single  Mode 

9 

- 

- 

- 

- 

5.6 

9.5 

- 

- 

Multi  Mode 

-32.6 

-188.2 

5.9 

4.  1 

5.4 

15.  1 

-0.5 

•24.4 

Ortho  Moments 

-32.6 

-186.  1 

6.2 

4.  1 

5.4 

15.3 

-0.5 

-24.4 

Air  Load 

-0.8 

0.4 

-0.4 

3.2 

-1.0 

2.3 

-0.7 

0 

Total  Deflection 

-35.9 

221.8 

-6.3 

158.9 

-3.7 

72.  1 

-4.2 

-72.3 

Singl e  Mode 

10 

- 

- 

- 

- 

- 

- 

7.4 

0.4 

Multi  Hode 

-256.4 

-391.5 

-9.  4 

24.3 

6.7 

30.9 

7.6 

-11.3 

Ortho  Moments 

-255.6 

-335.  1 

-9.4 

21.9 

6.4 

31.6 

7.6 

-11.8 

A  comparison  of  generalized  air  leads  deduced  from 
measured  lift  distributions  with  generalized  air  loads  obtained  from 
measured  moment  distributions  should  be  significant,  provided  the 
differences  are  not  lost  in  the  measurement  errors.  Unfortunately, 
this  appears  to  be  the  case  with  the  available  data,  and  a  definitive 
comparison  of  the  various  moment  reduction  techniques  cannot  be 
made.  This  is  demonstrated  by  a  plausible,  albeit  imprecise, 
analysis  of  the  air  load  and  moment  data  presented  in  the  next 
section. 
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4.  ACCURACY  REQUIREMENTS  ON  MEASUREMENT 
OF  BENDING  MOMENTS 


*» 


The  numerical  results  presented  in  the  previous  section  failed 
to  indicate  consistent  differences  among  the  results  that  would  permit  the 
definitive  evaluation  of  the  various  methods  for  analyzing  the  bending - 
moment  measurement  data.  It  was  postulated  that  these  differences 
might  have  been  lost  in  the  inaccuracies  of  the  basic  measurements  of 
both  the  air  loads  and  the  moments.  The  mannei  in  which  this 
possibility  was  investigated  is  described  in  this  section  and  proceeded 
as  follows: 

1.  An  estimate  was  made  of  the  threshold  and  accuracy 
of  the  measured  lift  distributions  presented  in  part  (b) 
of  Table  II  and  of  the  generalized  air  loads  derived  from 
them. 

2.  The  theoretical  maximum  bending  moment  per  unit 
generalized  air  load  in  each  mode  and  at  each 
frequency  was  determined.  These  quantities  multiplied 
by  the  appropriate  values  obtained  in  Step  1  constituted 
the  errors  in  the  bending  moments  that  would  be 
consistent  with  the  estimated  errors  in  the  measured 
lift  distributions. 

3.  Th®  requirements  derived  in  Step  2  on  the  threshold 
and  accuracy  of  bending-moment  measurements  were 
compared  with  estimates  of  the  corresponding  quantities 
in  the  actual  moment  data  as  presented  in  part  (a)  of 
Table  II. 

4.  On  the  basis  of  Step  3,  an  assessment  was  made  of  the 
validity  of  an  evaluation  of  analysis  techniques  which 
depends  upon  the  available  data. 
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THRESHOLD  AND  ACCURACY  OF  GENERALIZED  AIR  LOADS 


CALCULATED  FROM  MEASURED  LIFT  DISTRIBUTIONS 

The  lift  distributions  presented  in  References  12  and  14  were 
obtained  by  integrating  measured  pressure  distributions.  On  the  basis 
of  the  discussions  in  these  references,  the  accuracy  of  the  lift 
distributions  can  be  expected  to  be,  at  best,  within  +  2%  of  full  scale. 

It  is  (optimistically)  assumed  that  the  probable  error  is  constant 
(independent  of  radial  and  azimuthal  position)  and  equal  to  +  0.  2  pound 
per  inch.  This  value  is  only  about  2%  of  maximum  measured  lift  at 
the  inboard  radial  station.  Finally  these  estimates  must  be  interpreted 
in  terms  of  probable  errors  in  harmonic  content  at  each  radial  station. 
It  can  be  shown  that  the  corresponding  error  and  threshold  for  each 
harmonic  component  of  the  lift  are  both  equal  to  +  0.  1  pound  per  inch. 
The  scale  of  Figure  1  is  such  that  ..he  abscissa  lines  have  thicknesses 
corresponding  to  the  threshold  and  the  circles  drawn  about  the  data 
points  have  radii  equivalent  to  about  0.  1  pound  per  inch.  It  is  evident 
from  Figure  1  that  the  air-load  distributions  at  harmonics  beyond  the 
fifth  are  not  well  defined  for  this  case  [see  part  (c)  of  Table  II  ]. 

Estimation  of  the  accuracy  and  the  threshold  values  of  the 
generalized  air  loads  requires  an  assumption  about  the  manner  in 
which  the  lift  distribution  errors  are  distributed  with  respect  to  the 
particular  mode  shapes  considered.  For  the  purposes  of  the  present 
discussion  it  is  assumed  that  the  lift  errors  are  distributed  in  such  a 
way  that  the  generalized  air-load  error  is  maximized.  This  conserva¬ 
tism  is  at  least  partly  compensated  for  by  the  optimistic  assumptions 
regarding  the  accuracy  of  the  lift  distributions. 

The  maximum  error  in  the  generalized  air  load  is  given  by 

A(GJ.)nJk  =  *  f  {6k\\<t>n\cLr 

,  ro 

where 

=  error  in  the  k,ih  harmonic  component  of  the  air  lead 
measurement  at  each  radial  station,  pounds  per  inch 
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value  of  the  normalized  mode-shape  deflection  for  the  n 
mode. 


Similarly,  the  maximum  threshold  generalized  air  load  is  given  by 


where 

Ak  -  threshold  value  of  the  jcih  harmonic  component  of  the  air 
load  measurement  at  each  radial  station,  pounds  per  inch 

The  values  obtained  for  these  quantities  were  virtually  independent 
of  mode  shape  and  were  found  to  be 

A  (G.  A.  )  2;  +  12  pounds 

S  (G.  A.  )  +  12  pounds 

based  on  the  previously  obtained  values 

<5^  +  0.  1  pound  per  inch 

^  +  0.  1  pound  per  inch. 

A  review  of  the  generalized  air  loads  presented  in  Table  III 
in  light  of  these  possible  errors  shows  that  many  of  the  tabulated  values 
are  probably  highly  inaccurate.  This  is  demonstrated  in  Table  VI, where 
all  those  values  given  in  Table  III  which  are  less  than  the  estimated 
threshold  have  been  deleted  as  having  no  meaning  (indicated  by  the 
shaded  boxes)  and  where  the  possible  percentage  errors  in  the  remain¬ 
ing  values  are  as  indicated.  Almost  all  the  generalized  air  loads  for 
the  harmonics  above  the  fourth  are  less  than  the  estimated  threshold, 
and  many  of  those  indicated  for  the  lower  harmonics  are  subject  to 
large  errors.  Table  VI  indicates  that  the  evaluations  of  generalized 
air  loads  from  the  lift  distributions  are  estimated  to  be  in  error  by 
less  than  100%  for  only  the  following  cases:  Mode  1 :  k  =1,  2,  3; 

Mode  2:  k  =  l,  2,  4;  Mode  3:  fc  =  1,  2,  4;  and  Mode  4:  k  =1,  2,  3. 


f  |A, ,\\<t>n\dr 
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TABLE  TL 

ESTIMATED  PERCENTAGE  ERRORS  IN  THE  GENERALIZED  AiR  LOADS 
EVALUATED  FROM  MEASURED  PRESSURE  DISTRIBUTIONS 


r=— 

Harmonic 

Order 

k  =  oD/n. 

Mode  1 

ayn=  2.72 

Mode  2 
cOj/0-z:  4.95 

Mode  3 

«ya=  7.84 

Mode  4 

^/a  =  13.  1 

B 

2* 

II* 

oo 

19* 

2% 

CD 

9* 

83* 

1 

9* 

i3* 

30* 

16* 

fj|§§jiS 

9* 

7% 

21* 

34* 

3 

28% 

HK59H 

62* 

■ 

\2% 

- 

35* 

B 

Cosine 

28% 

57* 

96* 

- 

- 

64* 

81* 

- 

5 

32* 

.. 

- 

B 

- 

- 

- 

- 

B 

Cosine 

- 

96* 

- 

- 

32% 

- 

55* 

- 

B 

82% 

- 

- 

- 

~ 

- 

- 

8< 

- 

- 

- 

- 

- 

- 

9 

mm. 

- 

- 

- 

- 

- 

- 

10 

- 

- 

- 

- 

- 

- 

- 

MAXIMUM  MODAL  BENDING  MOMENT  PER  UNIT  GENERALIZED 
AIR  LOAD 


The  tip  deflection  (sine  or  cosine  component)  for  a  corres¬ 
ponding  unit  generalized  air  load  is,  from  equation  (7), 

_  _ l _ 

~  (uJnz  -Jo2az)  (u!/g)n 

and  the  accompanying  modal  spanwise  moment  distribution  is 


Mn  fcM  =  fnk. 


Since  the  maximum  moment  occurs  at  the  station  where  the  modal 


moment  per  unit  tip  deflection  is  a  maximum,  ??7„(max.)  ,  the  above 

equations  give  the  following  expression  for  the  maximum  modal  bend- 
ing  moment,  Mnk  ,  for  a  unit  generalized  air  load: 

7nn  ( man) _ 


M. 


nk. 


(u>Z-k.zflz)(w/<3)n 


(26) 


Values  corresponding  to  equation  (26)  are  given  in  Table  VII. 


TABLE  VI I 

MAXIMUM  MODAL  MOMENTS  FOR  UNIT  GENERALIZED  AIR  LOADS 


Harmonic 

Order 

ix  _ 

K  “  n. 

Mnk 

for  Unit  Generalized  Air  Load, 

inch-pounds 

Mode  1 

o01 

7T  =  2-72 

Mode  2 

U>z 

t r  =  *-96 

Mode  3 

COj 

7T  =7-84 

Mode  4 

cOd 

~JT  ~  l3, 1 

1 

5.51 

5.7  i 

5.21 

2.89 

2 

10.39 

6.54 

5.49 

2.95 

3 

21.90 

8.65 

6.00 

3.04 

4 

4.09 

15.72 

6.93 

3.  17 

5 

2.00 

304.0 

8.64 

3.37 

6 

1.23 

1  1.75 

12.38 

3.64 

7 

0.85 

5.50 

25.34 

4.02 

8 

0.62 

3.41 

122.3 

4.58 

9 

0.1*8 

2.38 

16.09 

5.43 

10 

0.38 

1.28 

8.  16 

6.85 

BENDING-MOMENT  MEASUREMENT  ACCURACIES  CORRESPONDING 
TO  LIFT  DISTRIBUTION  ACCURACIES 


The  accuracy  and  threshold  values  required  fr^m  bending- 
moment  measurements  to  obtain  generalized  air  loads  to  the  same 
accuracy  as  those  obtained  from  the  lift  distributions  can  now  be 
estimated.  The  appropriate  equations  are 


LM(n.k)  =  Mnjc  [A  (&.  A.)n\ 

(27) 

§M(n,k-)~  Mnk  [SCO.A.)^  • 

(28) 
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The  threshold  value  is  calculated  whenever  the  generalized  air  load 
in  Table  III  is  less  than  S  (G.  A. ).  12  pounds.  Results  are  shown  in 
Table  VIII. 


TABLE  IKE 

BENDING-MOMENT  MEASUREMENT  ACCURACY 
CORRESPONDING  TO  LIFT  DISTRIBUTION  ACCURACY 


Harmonic 

Requ  i  red 

Bending  Moment  Accuracy,  inch-pounds  -f 

Order 

Mode  | 

Mode  2 

Mode  3 

Mode  4 

K~ 

co/a 

<.4/1 1  =2.72 

a)z/n  =4.96 

ocJ3/a-  7.84 

tc)v/n.  =  13.  1 

,  J 

Cosine 

+66 

+68 

+63 

+35 

1 

Sine 

+66 

+68 

+63 

+35 

2' 

Cosine 

+  125 

+78 

+66 

+35 

Sine 

+  125 

+78 

+66 

+35 

Cosine 

+263 

+  i  04 

+72 

+37 

O 

Sine 

+263 

+  104 

+72 

+37 

,  J 

Cosine 

+49 

+  189 

+83 

+38 

*1 

Sine 

+49 

+  !89 

+83 

oo 

CO 

♦H 

5  ■ 

Cosine 

+24 

+3650 

+  104 

+41 

Sine 

+24 

+3650 

+41 

J  Cosine 

+  15 

+  141 

+  149 

+44 

L  Sine 

+  15 

+  141 

+  149 

+44 

7 

Cosine 

+  10 

+66 

+304 

+48 

.Sine 

+  10 

+66 

+304 

+48 

8J 

Cosine 

*r 

+41 

+  1470 

+55 

.Sine 

+7 

+41 

+  1470 

+55 

9 

Cosine 

+6 

+28 

+  193 

+65 

.Sine 

+6 

+28 

+  193 

+65 

i  Cosine 

+5 

+  15 

+98 

+82 

'  J 

[Sine 

+5 

■M 

+  15 

.... 

+98 

+82 

Table  VIII  indicates  the  accuracy  ana  threshold  requirements 
on  the  measurement  of  bending  moments  if  the  generalized  air  loads 
are  expected  to  be  determined  with  an  accuracy  of  ±  12  pounds  and  a 
threshold  of  ±12  pounds.  The  numbers  enclosed  in  shaded  boxes 
correspond  to  values  within  the  threshold.  Note,  again,  that  the 
moment  values  apply  at  the  station  at  which  the  particular  mode  has 
its  maximum  value  and  that  greater  accuracy  would  be  required  at 
all  other  stations. 
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At  this  point  it  should  be  recalled  that  the  preceding  arguments 
in  this  section  were  biased  to  reduce  the  requirements  on  bending- 
moment  accuracy.  The  maximum  estimated  error  was  used  with 
respect  to  obtaining  generalized  air  loads  from  the  pressure  distribution 
measurements,  and  the  maximum  mode  moment  was  used  in  estimating 
moment  measurement  requirements. 

THRESHOLD  AND  ACCURACY  OF  THE  BENDING-MOMENT 
MEASUREMENTS 


As  indicated  in  the  data  presented  in  part  (a)  of  Table  II,  the 
oscillatory  bending  moments  for  the  flight  condition  investigated  ranged 
from  -336  inch-pounds  3495  inch-pounds  at  —  -  0.  15  to  -6942 
<M<  5097  at  =  0.  65.  A  range  of  ±4000  inch-pounds  is  taken  for 
purposes  of  this  discussion.  The  accuracy  is  estimated  to  be  3%  and 
the  threshold  approximately  1%  of  the  reference  value.  These  give 
an  estimated  accuracy  of  ±120  inch-pounds  and  a  threshold  of  ±40 
inch-pounds. 

ACCURACY  OF  GENERALIZED  AIR  LOADS  OBTAINED  FROM 
BENDING-MOMENT  MEASUREMENTS 


A  comparison  of  the  bending-moment  measurement  accuracies 
estimated  above  with  the  requirements  presented  in  Table  VIII  shows 
that  the  generalized  air  loads  could  be  derived  from  the  bending 
moments  with  an  accuracy  at  least  as  good  as  from  the  lift  distribution 
for  only  the  following  points: 

Mode  1 :  k  =  2,  3 
Mode  2;  k  =  3,  4,  5,  6 
Mode  3:  k  =  5,  6,  7,  8,,  9 
Mode  4:  None 

Only  at  these  points  can  the  bending-moment  data  be 
considered  as  good  as  the  lift  distribution.  An  evaluation  of  the  various 
methods  of  analyzing  the  bending  moment  by  comparison  with  results 
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based  on  the  lift  distribution  is,  therefore,  possible,  at  most,  at  the 
above  points.  It  was  shown  previously,  however,  that  the  generalized 
air  loads  can  be  derived  from  the  lift  distributions  with  reasonable 
accuracy  at  only  certain  of  these  points.  Consequently,  a  valid 
comparison  is  possible  at  only  the  following  conditions: 


Mode 

1: 

k  =  2,  3 

Mode 

^  • 

■<* 

n 

Mode 

3: 

None 

Mode 

4: 

None 

Even  of  these,  the  expected  accuracies  of  the  reference  generalized 
air  loads  for  the  first  mode,  third  harmonic  and  the  second  mode, 
fourth  harmonic  are  not  very  good  (see  Table  VI).  A  definitive 
judgment  of  the  analysis  methods  cannot  be  made  on  the  basis  of  the 
extremely  small  sample  that  remains.  The  results  presented  in 
Tables  III  and  IV  are,  therefore,  not  indicative  of  the  relative  merits 
of  the  various  analysis  techniques  since  the  manner  in  which  the  data 
were  reduced  v/as  not  based  on  considerations  of  the  validity  of  the 
measurements. 

The  whole  discussion  of  accuracy  was  predicated  on  the 
assumption  that  the  blade  natural  frequencies  were  known  without 
error.  This  assumption  could  not  be  tested  in  the  present  program. 
Note  that  the  ability  to  detect  higher  harmonic  air  loads  from  bending 
moments  depends  strongly  on  the  amplification  of  the  bending  response 
near  resonance.  Errors  in  natural  frequency  determination  could, 
therefore,  degrade  the  results. 

Finally,  it  is  apparent  that  at  least  an  order  of  magnitude 
improvement  in  bending-moment  measurements  is  necessary  if 
approximate  generalized  air  loads  are  to  be  obtained  up  to  the  eighth 
harmonic.  It  is  probable  that  such  an  improvement  cannot  be  obtained 
without  the  development  of  specialized  filtering  and  recording 
equipment. 
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5. 


GENERALIZED  AIR  LOADS,  GENERALIZED  FORCES, 
AND  AIR- LOAD  DISTRIBUTIONS 


The  analysis  techniques  discussed  thus  far  in  this  report 
have  been  concerned  with  the  reduction  of  either  air-load  measure¬ 
ments  or  bending -moment  measurements  to  modal  tip  deflections  or, 
what  is  equivalent,  to  "generalized  air  loads.  "  Without  the 
introduction  of  any  assumptions  relative  to  the  air -load  distribution, 
the  generalized  air  load  constitutes  the  only  possible  interface  for 
comparing  results  based  on  the  two  different  measurements.  Recall, 
however,  that  the  generalized  air  load  is  the  integral  over  the  span  of 
the  product  of  the  mode  shape  times  the  lift  distribution  at  the 
harmonic  and  includes  the  lift  components  due  to  be  iding  slope 
and  bending  velocity.  As  a  consequence  of  this  dependence  on  the 
blade's  elastic  deformation,  the  generalized  air  load  does  not  provide 
the  blade  designer  with  the  basic  information  which  he  can  extrapolate 
to  the  analysis  of  other,  even  geometrically  similar,  blades.  Neither 
the  air-load  measurements  nor  the  bending -moment  measurements  can, 
by  themselves,  fulfill  this  need,  regardless  of  the  analysis  procedures. 
Therefore,  the  concept  of  the  "generalized  force"  (References  1,  2, 
and  5),  in  which  the  damping,  intermodal,  and  interharmonic  aero¬ 
dynamic  force  coefficients  are  assumed  known,  was  introduced.  It 
was  postulated  that  the  generalized  forces  for  one  rotor  blade  together 
with  the  assumed  motion-dependent  loads  might  then  be  used  to 
estimate  the  effects  of  the  air  loads  on  other  rotor  blades,  provided 
they  are  structurally  and  geometrically  similar.  The  relationship 
between  generalized  air  load  and  generalized  force  can  be  seen  in 
the  following  development. 

In  the  terminology  of  equation  (1),  the  spanwise  distribution 
of  aerodynamic  loading  is  approximated  by 

~  c(r)(nr+Vsinnt)(^-  +  V^LcosPd')  (29) 


where  t^Cr.t)  is  the  aerodynamic  loading  that  depends  only  on  radius, 
r  ,  and  time,  t  ,  (i.  e. ,  independent  of  blade  response),  and  the  last 
term  is  the  aerodynamic  component  that  depends  on  blade  bending 
velocity  and  bending  slope.  As  in  Section  1,  the  bending  deflection  is 
represented  by 


OO  OO 


«=/ 


and  the  aerodynamic  loading  is  given  by 

=  Z 

Again,  the  principle  of  virtual  work  is  invoked  and  the  orthogonality 
condition  on  the  natural  bending  mode  shapes  is  used  so  that,  finally, 
the  equation  expressing  the  equilibrium  of  transverse  forces  takes 
the  form  derived  in  Reference  5: 


(tajk  coshCit  -h  ia  sin  kCli  )  . 

'  C  5 


(31) 


E  (<W  - 

Ic.  —  O 


CO  -  , 

Ty2  fasC  k- 2)lli  -  cos(k*2)  fli) 

+£lV {cHfnn  ynkc  (cos(k+i) n-t  4-  cos(k-f)m ) 

-  JTlVVj  ,CMwn  $nkc  (cos(k-f)  fit  -  cas(k+  l)jTit  j 


(32) 


+  q,nk  cosh  fit 

~  mn  r  ”  5 


OO 


=  Z  Gmkr  cos  kSht 

k  =  0  c 


and 
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J  V\CHmn  $„ks  (si n  (kt2)nt  -  sin  (k -2) at) 
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,CM 
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-/  rn(r)  (p^(r)a'r 

uo 

=  f*c(r)  <pjr)  *ig±  dr 
Jo 

=  J'  rc(r)  <pm( r)  -  dr 

I*  c(r)  <pm(r)  <pn(r)dr 
<3 

f  rc(r)  <pm(r)(f>n(r)dr 


^mk. 


7mkj  c 


=-  f  4  kJr)  <f>m(r)  dr 
fR’tal  (r)<Pm(r)dr 


(33) 


(34) 
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The  generalized  force,  G  ^  ,  is  thus  seen  to  be  the  integral  over  the 
span  of  the  product  of  the  mih  mode  shape  times  the  k.th  harmonic  of 
only  that  part  of  the  aerodynamic  loading  which  is  independent  of  the 
blade  deformation.  The  generalized  force  does  not,  therefore, 
represent  the  effect  of  the  total  load  acting  on  the  blade,  and  the 
estimate  of  the  motion-dependent  air  loads  must  be  introduced  in  both 
the  evaluation  of  the  generalized  force  from  measured  data  and  its 
application  to  another  blade.  The  value  of  the  concept  of  the  generalized 
force  is  limited,  therefore,  by  the  error  introduced  in  the  approximation 
of  that  part  of  the  air -load  distribution. 

In  equations  (32)  and  (33)  for  the  cosine  and  sine  components 
of  the  generalized  force,  those  terms  involving  and  CM-mn  relate 

to  intermodal  coupling,  and  terms  involving  ic  ±1  or  fc  ±2  are  the 
interharmonic.  coupling  effects.  Even  in  this  relatively  complete 
representation,  the  unsteady  aerodynamic  effects  have  been 
approximated  by  quasi-steady  aerodynamics,  and  the  effects  of  reversed 
flow  have  been  neglected.  The  intermodal  and  interharmonic  aero¬ 
dynamic  coupling  effects  are  neglected  completely  when  the  generalized 
force  derives  from  the  Single-Mode  method  of  analyzing  bending- 
moment  measurement  data. 

The  discussion  above  is  in  terms  of  the  history  of  the 
generalized  force  idea..  It  was,  in  its  time,  a  useful  approach.  It  is 
reemphasized  here,  however,  that  the  errors  introduced  through  the 
assumptions  with  respect  to  the  distribution,  amplitude, and  phasing  of 
the  motion-dependent  forces  have  not  been  evaluated  and, hence,  the 
whole  scheme  must  be  used  with  caution.  An  insidious  situation  can 
develop.  For  example,  generalized  forces  derived  from  data  on  one 
rotor  and  applied  to  another  "similar"  rotor  could  give  perfect  results 
(within  the  accuracy  of  the  measurements)  if  the  two  were,  in  fact, 
aerodynamically  and  structurally  identical  and  if  tne  same  assumed 
motion-dependent  loads  are  used  in  both  cases.  But  this  same  result 
would  have  been  achieved  by  any  arbitrary  division  of  the  loads 
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between  motion-dependent  and  motion-independent  components  as  long 
as  the  same  division  was  observed  in  the  reconstruction.  The  limits 
of  applicability  of  schemes  for  deducing  generalized  forces  from 
generalized  air  loads  are  not  known,  and  experimental  data  to  establish 
the  required  "similarity”  are  not  available. 
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6.  CONCLUDING  REMARKS 


Theoretically,  no  other  set  of  basic  functions  can,  with  the 
same  finite  number  of  terms,  fit  the  bending-moment 
measurements  better  than  the  orthonormal-moment 
distributions  derived  in  this  report. 

"Generalized  air  loads"  constitute  the  only  interface  for 
comparing  analytical  procedures  based  on  bending -moment 
measurements  and  air -load  distribution  measurements. 

Accuracy  a.nd  threshold  limits  on  available  bending -moment 
and  pressure  measurements  preclude  the  definitive  evaluation 
of  the  various  methods  of  analyzing  moment  data  and, 
therefore,  the  practical  verification  of  Item  1,  above. 

Neither  the  bending -moment  measurements  nor  the 
air-load  distribution  measurements  are  sufficient  in 
themselves  to  provide  (without  assumptions  relative  to  the 
motion-dependent  aerodynamics)  a  semiempirical  method  for 
estimating  blade  exciting  forces  for  design  purposes. 

Whether  bending -moment  measurements  or  air- load  distri¬ 
bution  measurements  comprise  the  basis  for  the  evaluation 
of  the  generalized  air  loads,  additional  knowledge  of  the 
deformation-dependent  part  of  the  air  loads  is  required  in 
order  to  then  derive  the  generalized  forces. 
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7.  RECOMMENDATIONS 


1.  It  is  recommended  that  a  method  be  developed  for  evaluation 
of  the  motion-dependent  aerodynamic  forces  by  using  all  the 
experimental  information  available  (lift  distributions,  blade 
root  motions,  and  moment  distributions). 

All  the  following  recommendations  are  contingent  on  the  success  of 

Item  1 . 

2.  An  estimate  of  the  relative  importance  of  sources  of  errors 
should  be  obtained.  To  this  end,  it  is  recommended  that  a 
numerical  experiment  be  conducted.  Blade  structural 
characteristics  and  air  loads  could  be  postulated  and  used  to 
determine  moment  distributions,  natural  frequencies,  root 
slopes,  etc.  These  generated  data  then  could  be  treated  as 
errorless  measurements  from  which  air-load  distributions 
could  be  deduced.  Comparison  of  these  deduced  loads  with 
those  originally  specified  would  furnish  a  measure  of  the 
accuracies  obtained  with  different  reduction  techniques  (i.  e.  , 
bases),  dimensions  (number  of  modes),  etc.  Further,  the 
generated  data  could  be  perturbed  with  known  errors  and  an 
estimate  could  be  made  of  the  manner  in  which  these 
propagate.  An  effort  should  be  made  to  interpret  the 
findings  of  the  error  study  in  terms  of  (1)  the  accuracy 
with  which  natural  frequencies,  mode  shapes,  and  bending 
moments  must  be  known  and  (2)  preferred  gage  locations. 

3.  Development  of  specialized  filtering  equipment  is  recommended 
if  it  is  the  decision  of  USAAVLABS  to  utilize  cne  (or  more)  of 

the  reduction  techniques  discussed  in  this  report.  The  ) 

purpose  of  this  equipment  would  be  the  improvement  of  , 

measurement  accuracy.  This  could  be  done,  for  example, 
by  the  introduction  of  properly  phased  and  accurately 
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calibrated  signals  to  extract  most  of  the  zeroth,  first,  and 
second  harmonic  components  from  the  gage  signals.  The 
resulting  net  signals  then  could  be  amplified  to  display 
higher  harmonic  signal  contents  with  standard  recording 
equipment.  Use  of  the  device  in  model  tests  and  prototype 
flight  tests  would  enhance  the  structural  data  obtained. 

An  exploratory  study  of  the  use  of  edgewise  strain  measure¬ 
ments  to  construct  drag-load  distributions  is  recommended. 
The  initial  effort  should  be  devoted  to  obtaining  estimates  of 
measurement  accuracy  requirements. 

An  exploratory  study  of  the  feasibility  of  including  flapwise- 
edgewise-torsional  coupling  effects  is  recommended.  Such 
effects  might  not  be  important  with  respect  to  interpretation 
of  helicopter  rotor  structural  moment  distributions,  but  they 
could  be  important  to  the  interpretation  of  data  obtained  on 
highly  twisted  propellers. 
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